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I. INTRODUCTION 



Nowadays there is such an original situation in the theory of self-organizing systems. On the one hand, the synergetic 
concept has been developing successfully for more than twenty years. It allows us to explain the self-organization 
(ordering) of open system subjected to the environment disorder [jD. On the other hand, the phenomena of the 
self-organized criticality such as an avalanche motion of the sand grains on inclined surface (sandpile model) E|, 
intermittency in biological evolution H, earthquakes and forest fire, pinning in the random medium etc. (see |^) 
have been actively investigated for about ten years. However, in spite of the fact that the mission of both synergetics 
and theory of self-organized criticality is to explain the same phenomenon - the self-organization, they develop 
independently. It became so as within the framework of the synergetic approach the single statistical ensemble 
(formation of an avalanche) is investigated, whereas the examination of self-organized criticality models Q]- ^] is 
reduced usually to study of evolution of hierarchical avalanche ensemble [|J and is based on numerical methods and 
scaling representations || . 

Being the object of synergetics the formation of a single avalanche still remains an open question. It is evident 
that its reviewing is the only way to set connection between specified directions. Since phenomena of self-organized 
criticality are caused by avalanche ensemble behavior, and certain closed region of the state space corresponds to 
each avalanche, the standard formulation of the synergetic problem requires a nontrivial extension - while ordering 
self-organizing system we have to describe not only a symmetry breaking, but the ergodicity breaking that induces 
the clusterisation of phase space as well. In our paper a solution of this problem is suggested. In so doing we will 
investigate only the self-organization picture digressing from avalanche ensemble consideration (see 

The first attempt of field description of an avalanche formation was undertaken within the framework of the 
one-parameter approach based on nonlinear diffusion equation - thus, the feedback between open system and 
environment was not taken into account. Recently, two-parameter models were suggested, where the environment is 
represented either by control parameter ||, or conjugate field j[o|. According to || avalanche- like subcritical mode is 
formed in case of adiabatic relationship of characteristic times and first-order transition mechanism. The approach of 
Ref. allowed us to determine the critical indexes representing a scaling behavior of a self-organizing system within 
the framework of mean field approximation. Below we propose the generalized self-consistent scheme that is taken 
into account total number of the freedom degrees. Due to this we obtain not only mutually supplementary results 
P,|lO[ but also the complete analytical description of an avalanche formation. 

The suggested approach is based on synergetic generalization of the thermodynamic theory of phase transitions. 
The main feature of this theory is that within the conserved system (thermostat) a subsystem is segregated to 
represent a hydrodynamical mode whose amplitude is qualified as the order parameter pi[ that determines the 
state of the subsystem. Here it is assumed that thermostat influences upon value of the order parameter rj both 
thermally - by varying of a control parameter S, and immediately - by varying the field h conjugated to the order 
parameter (in case of magnet the values rj, h, S mean the magnetization, the magnetic field and the entropy). The 
distinctive peculiarity of the thermodynamic approach is that a postulation of the one-sided influence of thermostat 
on the ordering subsystem is declared, but not the reverse - the order parameter rj variation does not influence upon 
thermostat state parameters h and S. The synergetic approach considers the connection between an open subsystem 
and thermostat to be two-sided, so that the control parameter S and conjugate field h turn out to be the functions of 
order parameter rj. This connection manifests itself especially while the phase transition is described kineticaly. So, 
the standard picture meeting the Landau-Khalatnikov dissipative dynamics is realized in adiabatic approximation, 
when the relaxation time of the order parameter is much longer that corresponding times for conjugate field and 
control parameter . To this end the making use of Lorentz system, primarily suggested in order to describe the 
turbulent airflows | 13j| , happens to be rather convenient. 

Our paper is organized as follows. In sections || (basic expressions), III (Lagrange formalism), and VII (Appendix) 
is shown that the Lorentz system corresponds to the simplest Lagrangian of the supersymmetric field with components 
giving an order parameter rj, a conjugate field h, and entropy S. It is rather important that in this case a combination 
of the Grassmannian components of superfield plays the role of control parameter S in contradistinction to the 
usual field theory of a stochastic system fll3| ], where they are an auxiliary variables which have no physical meaning. 
The fact that the variables rj, h, S meet the vector of supersymmetric space is a reflection of the self-consistent 
behavior of the synergetic system (in contrast to the statistical field scheme Jul, where the superfield is only a 



convenient technical method). The study of superfield correlators is carried out in section IV. As is known the 
components of such correlators are not independent in ergodic state - the presence of supersymmetry causes the 
fluctuation-dissipation theorem which connects stated components H]. If the quenched disorder is appeared, an 
ergodicity loosing happens that breaks the supersymmetry in its turn. This leads to appearance of singular additives 
to correlators that define memory and nonergodicity parameters q, A. The basic result of our work is the defining 
quantities q and A dependencies on thermal and quenched disorders intensities. This allows us to find the conditions 
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of avalanches formation in the self-organized criticality mode appearance. The generic example of such a process is 
considered in section ^ as the flow of sand on inclined surface. It is shown that this process can be represented by 
the Lorentz system, considering the horizontal and vertical components of the sand grain velocities to be the order 
parameter rj and conjugate field h, and the tangent of a surface inclination angle as the control parameter S. Final 
section |^| is devoted to discussion of obtained results. It is shown that a parameter which determines the transition 
to time irreversible regime is given by ratio of a time of the quantum fluctuation to the one of order parameter. A 
critical value of effective interaction that bounds the domain of ordered state is found. 

II. LORENTZ SYSTEM 

To expound the microscopic scheme of the self-organization description let us study firstly the system consisting 
of Bozon and Fermion gases whose interaction is characterized by potential v. Within the framework of secondary 
quantization Bozons are described by the and 6; operators, satisfying the usual commutation relation: [bi, 6+] = <5; m , 
where I, m are the site numbers. The two- level Fermion subsystem is represented by operators o^, ai a , a = 1, 2, for 
which the anti-commutation relation {oia^a^p} = ^imS a p is fulfilled. The occupation numbers b^b k determine the 
Bozon distribution within k- representation that corresponds to the Fourier transform over lattice sites I. To represent 
the Fermi subsystem we should introduce the operator di = a^a^ determining the polarization with respect to the 
saturation over levels a = 1, 2, as an addition to the occupation numbers ni a = a^a; Q . As a result, behavior of the 
system under consideration is defined by Dicke Hamiltonian 



H = Y1 |(^i n ki + E 2 n k2 ) + u> k b+b k + ^v(b+d k - d+6 k ) J , (1) 



where the k-representation is used, Ex 2 are the Fermi levels energies, uj k represents the Bozon dispersion law and the 
imaginary unit before the interaction term v reflects the Hermitity property, the Planck constant is % = 1. 
The Heisenberg equations of motion corresponding to Hamiltonian p) have the form 

i k = -i^k&k + (v/2)d k , (2a) 



c?k = -iArfk + (w/2)6 k (^k2 - Tiki), (2b) 
Tiki = (v/2)(b+d k + d+b k ), (2c) 

n k2 = -(v/2)(b+d k + d+b k ), (2d) 

where the dot stands for a derivative with respect to time and the quantity A = E 2 — E\ is introduced. In case 
of resonance suggested the first terms in the right-hand sides of equations (pa|), ( ph| ) that contains frequencies Wk 
and A may be eliminated by extracting the multipliers exp(— iw k i) and exp(— iAt) in the time dependencies b k (t), 
d k (t), respectively. On the other hand, if the dissipation is taken into account, these frequencies obtain imaginary 
additions — i/r, ; , — i/r^ characterized by relaxation times r r) , t/j (here the conditions Im Wk < 0, Im A < reflect the 
'causality" principle). As a result, equations (|2a|), (2b) get the dissipative terms —b^/r^, —d^/rh, where t v is the 



relaxation time of Bozon distribution and is the Fermion polarization time. One can suppose that the dissipation 
influences onto the Fermi levels occupancies n ka (t) also. However, since the stationary values n ka ^ (and in case of 



external pumping n k2 > n kl ) the dissipative terms in Eqs. (pq), (2d) have much complicated form: — {n ka — n ka )/Ts, 
a = 1, 2, where is the relaxation time of the Fermion distribution over level. 
Now, let us introduce the macroscopic quantities: 

77k = (6+) = (6 k ), h k ee (d k ) = (d+), 

Sk = (n k 2 - n k i), 5k = (n k2 - n kl ), (3) 
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where the angular brackets mean thermodynamic averaging. Neglecting the correlation in distribution of particles 
over quantum states the Heisenberg equations (2), being contemplated by dissipative terms, result in the Lorentz 
system 

Tr,V = -r) + A v h, (4a) 
r h h = -h + A h r)S, (4b) 

t s S = (S° -S)- A sV h. (4c) 

Here in terms of one-mode approximation the dependence on the wave vector k is omitted and the constants defined 
by relationships 2A V = vt v , 2 Ah = VTh, As = 2vts are introduced. Equations (4) contain following seven constants 
- pumping parameter 5°, three relaxation times r and three coupling constants A. But as there are above relations 
between these last ones caused by interaction parameter v > 0, only five of them are independent. Since four of these 
fix the scales for quantities rj, h, S, i, so only the parameter of thermal disorder S° plays a substantial role and its 
value determines only the system behavior EL^2) . 

To analyze equations (4) we introduce the scales rj m , h mi S c defining ranges of the variation for the order parameter 
77, the conjugate field h and the control parameter S: 

i lm 2 = A h A s = T h T S v 2 , h- 1 = A n / Vm = (T v /2)(T h T S ) 1/2 v 2 ; 

S- l =A ri A h = 2- 2 T ri T h v 2 . (5) 
Then, using magnitudes 77, h, S normalized by 7] m , h m , S c values, we results Eqs. (4) in the form: 

t v t] = —77 + h, (6a) 

Th.h = —h + r]S, (6b) 

t s S = (S° - S) - rjh. (6c) 



In terms of adiabatic approximation r^ , t$ <C t j( the left-hand sides of equations (|6b[), ( J6q ) may be set equal to zero. 
Thus, we derive to a result 

h = S%/{\ + r l 2 ) 1 S^S'/il + v 2 ). (7) 

With the order parameter growth in physical domain 77 S [0, 1] the conjugate field increases and the control parameter 
decreases monotonically. If rj > 1 the h{rj) dependence is of decreasing shape, that corresponds to unstable state. 
Inserting ((?]) into ( |6a| ) we obtain the Landau-Khalatnikov equation 

T n f, - -dV/dr,, V^ 1 - {if - ^° ln(l + ry 2 )) . (8) 

Its form is defined by synergetic potential V(ri) having its minimum at point rjQ=(S°— l) 1 / 2 . Hence it is seen that 
stationary value of the order parameter rjQ ^ is realized under the thermal disorder conditions S° > 1 (5° > S c 
in usual units). Thus, the magnitude S c defined in last equality (ra) is the critical value of the control parameter. 
According to Eq. ^ in the stationary state one has ho — {S a — l) 1 ^, Sq — 1. The last equality implies that despite 
supercritical value S° > 1 of thermal disorder the system relaxes so, that the stationary value of the control parameter 
So = 1 is reduced to the critical one. 

The mentioned relaxation is provided by the negative feedback of order parameter rj and conjugate field h with 
control parameter S that is described by the last term of equation (pcf). This feedback, displaying the Le Chatelier 
principle for the self-organizing system, compensates the Sir]), /i(?7)thermostat's state parameters increase which 
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takes place, when this feedback is absent. On the other hand, the positive feedback of quantities r\ and S with h in 
( |6h| ) is the reason for the self-organization. It is obvious that the stationary state rjo, ho, So can be realized only under 
the condition of inverse subsystem influence, characterized by order parameter r], on the thermostat parameters h, S. 
It is worthwhile to note that inverting the signs of nonlinear terms in equations (6) causes the minus appearance in 
the right-hand side of the first equality ([?]), so that the susceptibility \ — dr//dh becomes negative and this case does 
not meet the stable state. 

Thus, the self-organization process takes place only if the negative feedback of order parameter rj and field h with 
control parameter S and the positive feedback of rj and S with h both exist. According to (6b), ( |6q ) such a choice 
of signs is determined by the fact that the negative feedback provides falling-down of the control parameter S in the 
course of time, whereas the positive one insures the field h growth. Further we shall show that the value S is reduced 
to the entropy and its decrease reflects the non-conservation of the self-organizing system for which the second law of 
thermodynamics is not fulfilled. The crucial role of the increasing character of field h is stipulated by the fact that 
the linear equation (6a) for the order parameter rj contains namely the field h. As a result, influence of the field h 
on the velocity of r)(t) increase, and on the self-organization process also, is direct, whereas the influence of control 
parameter is indirect. 

The described scheme of self-organization meets the second-order phase transition. To describe the first-order 
transition one ought to set the relaxation time t v of order parameter to be the function of its value rj ]l],[l2| . Such a 
scheme of an avalanche formation is represented in JtJ . 



III. LAGRANGE FORMALISM 



In the previous section we omitted, within the adiabatic approximation, the fluctuations of conjugate field h and 
control parameter S and, thus, we made it possible to reduce the Lorentz system (6) to the Landau-Khalatnikov 
equation (||). To form the Lagrange formalism one should accomplish the reverse transition supposing a fluctuation 
source £ appearance in (|^). If the nonhomogeneity is considered, the generic expression is reduced to the Langevine 



equation (see 15 for example) 

fj{v,t)-V 2 rj{v,t) =/(r,i) + C(r,i), (9) 

where r is the coordinate measured in units of the correlation lengths £ and t is the time related to the scale t v , 
the force / = — Vo'{rj), Vo = dVo/drj is defined by dependence Vo{rj) for the bare potential related to fluctuations 
intensity T. The term — V 2 ry in the left-hand side of Eq. (|9|) takes into account the spatial nonhomogeneity within 
the framework of Ginzburg-Landau model. The expression (m is valid for nonconserved order parameter, otherwise 
terms —V 2 rj, / obtain the additional operator —V 2 jlq ]. The fluctuational term is normalized by the white noise 
conditions 

(CM))=0, (C(r,t)C(r',0)=TO(r-r')^-t'), (10) 

which correspond to averaging over the Gaussian distribution with dispersion T. 

To construct the Lagrangian corresponding to the Langevine equation (|J), let us use the standard field scheme JTJ] 
based on a generating functional 




ZMr, t)} = (T\5(T,- V 2 rj - / - C) det 



(11) 



being the generalization of the partition function. Herein the continual product of <5-functions takes into account that 
the condition ([)]) ought to be satisfied with all values r, t and the determinant represents the Jacobian of ( to rj 
transition. We apply the Fourier transform for (5-function leading to appearance of the field f(r,t) and introduce the 
Grassmannian conjugate fields ip(r,t), i/5(r,i) into operation for integrated representation of the determinant. Thus, 
the equality ( |TT| ) assumes the canonical form 



z iv} = J p {V, 9, V'j ipjDipDtpDtp, Poce" s , S = J C(rj, </?»)drdt, (12) 
where the Lagrangian 

C = v{n - V\) + - VV) - <^ 2 /2 + <pV '(ri) + V a "(v)W> (13) 



5 



is measured in noise intensity T units. The form of the corresponding Euler equations is as follows: 

V - V 2 ?? = -V Q '(v) + <P, 



(14a) 



<p + V 2 <p = V "(Ti)<p + Vo m (v)W, (14b) 
i> - V a V - -Vb"(r#, (14c) 

$ + V 2 Tp = V "(rj)%j). (14d) 

The first one can be reduced to the Langevine equation (|^) by replacing the field tp by stochastic component Q. Since 
this equation corresponds to the maximum of a probability P distribution in ( |l2| ) , tp represents the amplitude of the 
most probable fluctuation of the field conjugated to the order parameter rj (the force / is the average value of this 
field). Obviously, the conditions {() = 0, tp =/= imply that during the self-organization the bare Gaussian distribution 
transforms from unimodal into bimodal form with its maximums at points ±tp. The signs distribution in front of the 
gradient terms in Eqs. (14) is quite remarkable: albeit the standard combination inherent in relaxation processes like 
diffusion is realized for components r], tp, but the fields tp, tp contain nonhomogeneity terms with opposite signs that 
means the autocatalitical increase of those components. It is shown in Appendix how the field equations (14) are 
reduced to the Lorentz system (6). 

The most elegant method to represent developed field scheme is to incorporate the components rj, tp, tp an d the 
generalized force <p = — 6Vo{r]} / Srj into supersymmetrical field 

$ = V + ip~X + XV> + XX0, (15) 

where Grassmannian coordinates \i X posses the same anti-commutation properties as the ip, Tp fields. To represent 
the Lagrangian ( |l3| ) supersymmetricaly one should first of all replace the bare potential Vo(r)) by renormalized one 
V(rj), presented in motion equation flA.7| ), ^ and get rid of gradient addends, committing transform to the variational 
derivatives V'{r]} = SV{rj}/5r] = dV(rj)/dr) — V 2 ?7, V{ rf\ = JV(rj)dr. Then, expressing fluctuation amplitude tp in 
terms of generalized force <p according to equality ( |A.9| ), we derive to the following form of the Lagrangian (p"3|): 

C = (V/2 + H>~ <t> 2 /^) + (-V'{ri}<p + i)V"{q}i) \ + V'Mv- (16) 

The last addend can be omitted as the total time derivative of ^{77} and in superfield representation (15) the 
Lagrangian ( |l6| ) assumes the canonical form: 

C = J A($)dxdx, A= {1/2)§DD$> + V{§). (17) 

Here the kinetic superenergy of the kernel A corresponds to the first bracket of expression ( |l6|) and the potential 



superenergy V($) - to the second one. The equalities ( A. 3 ), (A. 7), and (A.S) within the framework of the <l> 4 -model 
give 

V= 1 -^& + J-±^Z\ (18) 

where the anharmonicity parameter w > emerges because in contrast to the order parameter rj (see @) the superfield 
( |l5| ) cannot be scaled by the only magnitude r\ m . The supersymmetry group generators have the form 



It is shown in Appendix that the renormalization is caused by self-consistency of the superfield components (jl5|) . 
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The superequations of motion ensuing from the superaction extremum condition S{^(z)}— J A($(z))dz, z 
{r,t,x,x} read: 



(1/2)15, D]$ + V'{$} = 0, (1/2)[D, £>]$ = + ^x + XV- + XXV- 
Projecting Eqs. ( po|) onto basis vectors 1, x, Xj XX °f the superspace, we arrive at equations 

V = -V"{i 1 }cj) + V'"{r,}^, 



(20) 



(21a) 



<t> = -v'{vh 
i> = -v"{n}i>, 



(21b) 
(21c) 



i/j = V"{q}i). 



(21d) 



The last of this set can be obtained from (14c), (14d) by replacing Vq by V , and the re latio nship (211) gives the 
definition of the force 4>. Eq uatio n (21a) is obtained by time differentiating of equation (14a) and substituting the 



derivatives fj, (p from (A.E), (14b) into resultant expression. Thus, systems (14), (21) turn out to be equivalent with 
accuracy to the bare potential Vq(v) renormalization. However, though the equations of the first set are symmetrical 
concerning the time derivative order, in (21) this symmetr y is broken due to the transformation the fluctuation tp to 
the generalized force (j>. Comparing Lagrangians ( |13| ) and (|l6|) shows the fact that the above transformation provides 
the standard bilinear form of super lagrangian with respect to operators D and D. Easy to show that, if the gauge 
condition Z)$ — is satisfied the Grassmannian fields %p, ip are suppressed and the single combination XX replaces 
the pair of conjugate ones x, X- 111 this case the kinetic superenergy is linear relatively to the supergroup generator, 
moreover both pairs - ij, </> and ij, ip are allowable to be used as superfield components |fl7j| , However, as is shown 
in Appendix, at the description of the self-organization the behaviour of the Grassmannian fields represents entropy 
and, consequently, is essential. Therefore the gauge D$ = is not fulfilled, and one should prefer making use of the 
generalized force (f> to the fluctuation amplitude (p. 



IV. CORRELATION TECHNICS 



Let us introduce the supersymmetrical correlator 

C(z,z') = (*(*)$(/)>, z = {r,i,x,x}- (22) 

According to the motion equation ( ^o|) its bare component C^(z,z') meeting the potential = (1 — cr)$ 2 /2 
satisfies equality 

LUX)C ( ^(X,X')=S(X,X'), L=(l-a) + (l/2)[D,D], (23) 

where the transition to the time-spatial Fourier transforms is made and the supersymmetrical ^-function 5(x,x') = 
~{x ~ x')(x — x') is introduced. Taking into consideration definitions ( |l^ ) and expressions (l/4)[£>, D] 2 = —ui 2 , we 
obtain 

C (0) (x ^)- 1 ~^ x+ ^ /x/ ) + [( 1 ~^ +iu; ^ x/+ [( 1 ~ c^ )^' ia '^ /x ~ a ' 2 ^ /x, (24) 
' (1— a) 2 +LU 2 

where index uj and spatial dispersion were omitted for brevity. Equation (EJ) is an expansion in basis components 
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B o = -XX, Bi = -x X , 

T = l, T!=XXXV; (25) 
Fo = XX, Fx = xx' 

whose functional product satisfies the following multiplication rules: B 2 , = Bo, B 2 = Bi, Fq = Fo, F 2 = Fi, 
B Ti = Ti, BiT = T, TB = T, TTi = Bi, T1B1 = Ti, TiT = B , and the other multiplicands are equal to 
zero. Thus, Bq,i, T, Ti, Fo,i form the closed basis, and it is convenient to expand supercorrelator ( |22] ) in these 
components: 

C = ,g + B + 0-B1 +ST + sTi + G+F + GLFi. (26) 
Insertion ([l5]) into ( p2| ) derives to the coefficients 

9+ = 9- = -{ i n4>); 

s = (v 2 ), * = -<M 2 ); (27) 
G+ = -(M), G_ = (W), 



where the fact is considered that according to (A..9) the field (f> is exclusively imaginary. Thus, the magnitudes g± 
are reduced to advanced and retarded response functions of the order parameter rj to field <f> action; S and s are the 
autocorrelators of order parameter 77 and field cj), and the functions G± determine the correlation of Grassmannian 
conjugated fields "ip, ip. The statements 

3 £ 0) = = [(i - of + to 2 }- 1 , * (0) - -^ 2 [(i - °) 2 + u 2 ]-\ 

Gf = [(l-ffj+iw]" 1 , G (0) = [{l-a)-iuj]- 1 (28) 



are valid for the bare supercorrelator (|2~i|). Hence, according to ( A.9 ) the relation (rj<p)o=G_ is deduced that is a 
special case of the Ward identity Jl4j] . It means that the correlator of the Grassmannian fields is reduced to response 
function of the order parameter rj to the fluctuation ip. 

Expansion (^) allows us to treat supercorrelator ( p2| ) as the vector of the direct product of superspaces. Since 
the Fermi components Fo,i do not couple with the Bose ones Bo,i, T, Ti, we may use more compact basis, passing 
on from the field <f> to the fluctuation ip. To doing so, we will neglect the correlation of the self-consistent field 4> 
whose structure factor is s ~ uj 2 , when uj — > 0. Moreover, we will replace the response functions g± to field <p by 
corresponding functions G± for the fluctuation ip. Then the expansion ( p6|) takes such a compact form 

C = G+A + G_B + ST, (29) 

where the basis operators A = Bo + Fo, B = Bi + Fi are introduced that satisfy the multiplication rules: A 2 = A, 
B 2 = B, BT = T, TA = T (other multiplicands are equal to zero). 

In the issue the behavior of self-organizing system is described by the Lagrangian 

L = (ipri + tptp - ip 1 12) + {V'i^ip + 4,V"{r]}tl;). (30) 
We introduce the quenched disorder 



2 (#r) - 0) 2 - (A^y 



? = (A ; )2 (3D 



Its magnitude characterizes the field random scattering 4>(r) (the dash in (|3lj) stands for averaging over coordinate r, 
the (Aiy?) 2 = |<£ w= o| 2 is the square- mean fluctuation). If the quenched disorder is brought in the action, meeting the 
Lagrangian (pG), component squared in ip fluctuation takes the form 

1 f \^\ 2 ^-^ /5M|^| 2 dc. (32) 



2 J ' rw < 2vr 

Here we have neglected an integration over r and passed to Fourier transform over frequency uj. At the equilibrium 
disorder the field scatter <fi(r) is reduced to the square-mean fluctuation (Ai^) 2 so, that p — and expression j3^ ) 
has the canonical form —(1/2) J ip 2 dt. In the case of quenching one has p > and second addend in (|32|) leads 
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to renormalization of the bare supercorrelator ( pi] ) whose component gets the multiplier 1 + 2ttp 2 S(uj) in (|2 
Respectively, we find the operator L in the motion equation ( |23|) as follows: 

L = L + A + LB + LT; L± = (1 - a) ± iw, L = - [l + 2t7 P 2 5(lj)] . (33) 



To obtain equation that defines the supercorrelator ( |22|) one should multiply (|2C|) by $(z') and average the result 
over distribution -P{$} from (p^). Thus, we get the Dyson superequation 

C^L-S, (34) 

where within the framework of <E> 4 -model ( |l8| ) the self-energy superoperator S is defined by equality 

E(z,z')-(2/3)w 2 (l + 3a) 2 (C(z,z')) 3 , z = {r,t, X ,x}- (35) 

Herein the w > is the anharmonicity parameter related to temperature T and the condition J C(z, z)dz — is taken 
into account that follows from Eqs. ( p5| ) and (p9|). 

If the anharmonicity is omitted, E = 0, the components (|2g|), diverging at the point of transition into self- 
organization state (a — 1), are obtained from Eq. (f34|). Thus, the supersymmetrical field approach allows us to 
reproduce the main result following from the Lorentz system (6) by means of the linear approximation. In the general 
case the self-energy superfunction should be expanded similar the supercorrelator (p9|): 

S = S+A + £_B + ST. (36) 

Then, according to Eq. ( |33"| ) the superequation ( p4] ) is reduced to components 

G±~ 1 = [(1 - a) ±ilu] - S±, (37a) 

S=[l + 2ip 2 %) + E]G + G_. (37b) 

The explicit form of the expansion coefficients ( |36| ) is given by expression (|35|). In accordance with the su- 
percorrelators multiplication should be understood in usual meaning, not in functional one: (T(x,x')) 2 = ^(XjXO) 
Mx, X')T( X , x') = T( X , X')A( X , X') = Mx, X'), B( X , x')T(x, x') = T( X , x')B( X , x'\ = B (x,£), and the other multi- 
plicands are equal to zero. As a result, for the spatially homogeneous case we obtain from ( pq ) the following 

E±(t) = 2w 2 (l + 3<r) 2 S 2 (t)G±(t), (38a) 

E(t) = (2/3)w 2 (l + 3(i) 2 5 3 (i). (38b) 
At insertion of these expressions into the Dyson equation (37) we will be in need of their frequency representation 



that contains convolutions. To avoid such a difficulty let us use the fluctuation-dissipation theorem |14 1 



SV) = (2/u/)ImG±(u/), £(u/) = (2/o/)±mE ± (a/), (39) 

where cj' is the real frequency. Using the spectral representation of the complex frequency lo and integrating equalities 
©, we find 

5(t 0) = G±(w -> 0), E(t -> 0) = E±(w 0). (40) 
Since the G±(w — > 0) gives x i n the hydrodynamical limit — > 0, we have 

^^0) = X = G±(^^0), (41a) 

E±(w^0) = (2/3)u; 2 (l + 3a) 2 x 3 , (41b) 
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where the expression ( p8q ) is used in 
second order of the anharmonicity w. 

The equations (37), (38), and (41) describe the behavior of the self-organizing system completely. Particularly, they 
represent not only the ordering phenomena but the effects of ergodicity breaking and memory appearance as well. 
These effects manifest themselves in elongation of correlators 

G_(w) = -A + G_ (w), S(t)=q + S (t) (42) 

at the expanse of the Edwards- Anderson memory parameter q = (77(00)77(0)) and irreversible response A = \ — xo that 
is equal to difference of microscopic susceptibility \ = G-q(uj = 0) and macroscopic mag nitude Xo = G_(w = O).0> 

Now, let us insert the elongated correlators (E3) into expressions (38). Then the renormalized components of the 
self-energy superfunction read 

E±(t) = 2w 2 (l + 3a)V (-A + G±o(t)) + E± (t), 

E ±0 (t) = 2w 2 (l + 3a) 2 S (t)G ±0 (t)(2q + S (t)); (43a) 



41b). In contrast to |19] here the self-energy components E± contain only the 



£(rj) = (2/3)w 2 (l 



3a)V(<? 



35 (t)) + E (t), 



E (t) = (2/3)tt; 2 (1 + 3ct) 2 S 2 (t)(3g + 5 (t)). 



(43b) 



Here the no nlin ear terms with res pect to the correlators G±o, So are ga ther ed in addends E±0) So, in the second 
equality of (|43a|) , the relationship (41a) is considered, and in equations ( |43b| ) the addends containing Sq A ~ arc 
omitt ed. I f the memory is absent the first addends of E±(t), E(f) vanish. Inserting the Fourier transforms of equations 
(p2|), (43b) into the Dyson equation (37b), within w-representation we obtain the following relations: 



qo [1 - (2/3)7^(1 + iafxlq 2 ] = p 2 X l 



(44) 



(1 + E )G+G_ 



1-2t7j 2 (1 + 3a) 2 q 2 G+G. 



(45) 



The first one corresponds to the <5-like addend caused by the presence of memory, when u> = 0, and the second one 
meets the frequencies 10 — > 0. When ui = the characteristic combination G + G_ goes to G + G_ = \o an d the pole of 
the structure factor (EF 



2t« 2 (1 + 3 ( 7 c ) 2 (7 2 



(46) 



determines the point of ergodicity breaking c c , where Xq=X, Qo = Q- 

Now let us insert the Fourier transform of expression (43a) into equation (37a). With the help of equality (41b) we 
obtain the following equation 



GI 1 + 2w 2 (l + 3a) 2 q 2 G- = [(1 -a)- vu] - 2w 2 (l + iaf X 2 {q + x/3) 



(47) 



for the retarded Green function in the hydrodynamical limit to — > 0, where the relationship (41a) is taken into account. 
Hence the expression springs out for the microscopical susceptibility x=G_(w-^0) 



1 - (1 - a) X + (2/3)™ 2 (l + 3a) 2 x [(x + q) 3 - q 3 ] = 0. 



(48) 



The behavior of self-organizing system with quenched disorder is determined completely by equations (}44|), (|4€|), 
and (^8|). Herewith one should contradistinguish the macroscopical magnitudes qo, Xo and microscopical ones q, x 



' Let us point out the reverse sign of the irreversible response A as compared with the definition given for thermodynamical 
systems, where the ordering corresponds to low values of the noise intensity (temperature). 

3 ' The unified function G_(w) can be used to define susceptibilities \o, X considering that magnitudes xo=G_(a;=0) and 
X = G_ {to — » 0) correspond to equilibrium and out-of-equilibrium values. In so doing we should equip all correlators in Eqs. ( fl"c| ) 
and (41) with index and set uj — 0. 



10 



(the first correspond to frequency 10 = and the second to limit u> — ► 0). The distinctive feature of this hierarchy is 
that macroscopical values qo , xo depend exclusively on quenched disorder intensity p while the microscopical ones q, 
X - on thermal disorder a. Thereafter, to determine values qo, x.o the magnitude a should be considered to be equal 
to value <J c {p) on the line of ergodicity breaking, and the quenched disorder intensity takes the critical value p c {&) for 
the defining of q, x- Hence, equations (Bjj), ( [fq ) determine the macroscopical values qo, Xo and equation (^]) defines 
microscopical ones q and X- Herewith the external addition (f> ext to the self-consistent field <f> (hereinafter <j) ext = 
is supposed) , the dispersion p of field <p fixing the intensity of the quenched disorder (|Tj) , and the parameter of the 
thermal disorder a = S° / S c act as the state parameters. 

Combining the equalities (p4|), (fl6"|), we get the expression for the macroscopical memory parameter 

qo = (3/4) x /3(l + Sa c y V3 (p/w) 2/3 , (49) 

which increases with p growth. Inserting Eq. ( |49[ ) into Eq. (^6|) gives the macroscopical susceptibility 

Xo = 2 1 / 6 3- 1 / 3 «;- 1 / 3 (l + 3a c )-^ 3 p- 2 / 3 (50) 

that decreases with p growth. 

Fixing the memory parameter in equation ( [l8|) by expression (compare it with (|49|)) 

q = {3/4)^ 3 (l + 3<j)- 2 / 3 {p c /w) 2 / 3 , (51) 

we find microscopical values q, x dependence on a (Figs. 1, 2). The fact that, when the thermal disorder is small the 
functions q{u) and x( a ) are two- valued is distinctive for these dependencies. Since the susceptibility should increase 
near the ordering point the values q and Xj shown in Figs. 1, 2 by the dotted lines, are nonstable. There is an abruption 
at the ordering point a c , where dx/dc = oo. The growth of anharmonicity parameter w leads to the reduction of 
both susceptibilities x an d Xo- 

The ergodicity-breaking point is fixed by equation 

3,4 + (A+ l) 3 = 1 + 2(1 -a c )p- 2 , A = 2 5/6 3- 2/3 {l + 3a c ) 1/3 w 1/3 p-' l/3 , (52) 

obtained from condition x — Xo according to equalities (^8|)-(pT|). As is seen from Fig. 2 the quantity a c that meets to 
the ordering point is defined by the maximum value a c corresponding to the breaking of ergodicity. The dependencies 
cr c (p), ct c {p) which represent the phase diagram of the self-organizing system are shown in Fig. 3. The disordered 
state region corresponding to small values of a, and the position of it's boundary does not depend on p, whereas 
the ergodicity region is bound by small values of p. [j Comparing Fig. 3a with 3b shows that the behavior of self- 
organizing system is rather sensitive to the anharmonicity parameter w picked. The appropriate dependence a c (w) 
for the thermal disorder parameter that corresponds to ordering, is adduced in Fig. 4. The a c takes its maximum 
value a c = 1, when w = 0, and the growth of w causes the a c monotonous decrease until er c = at w c = 0.064. The 
maximum value p m = p c (a — 0) of the quenched disorder at the ergodicity breaking line changes with the value w 
according to dependence shown in Fig. 5. When w is small, the quantity p m oc w^ 1 ^ 2 increases infinitely and the 
ergodicity region disappears with growth of anharmonicity parameter above critical value w c . 

Nonergodicity parameter is defined by solution of Eqs. (^q)-(|5l|). At constant value of quenched disorder intensity 
p (Fig. 6a) the three regimes are possible. At small p macroscopic susceptibility xo exceeds microscopic quantity x( a ) 
for every values of a (see Fig. 2) , and system is always in nonergodic state. When parameter p reaches the values which 
exceeds threshold p c meeting the condition xo(p) = x( cr c), the nonergodicity parameter A = x — Xo takes the nonzero 
values only within the range a > a c . Microscopic susceptibility x for every values of a exceeds the macroscopic one 
Xo starting from value p m that meets the condition x.o(p) — x( a = 0) an d system is nonergodic always. At fixed 
value of thermal disorder a (Fig. 6b) the dependence A(p) is defined by infinite increase of macroscopic susceptibility 
Xo ~ p~ 2 / 3 in the region of weak quenched disorder p — > 0. As a result, nonergodicity parameter takes the nonzero 
values starting from critical value p c , and increases at further growth of quenched disorder monotonically. 

Analytic representation of A(a,p) dependence is possible only near the line of ergodicity breaking <J c (p). Setting in 
Eq. (48 ) < a — a c <§; a c , x — Xo + A, A <S x i n first order over small values (a — <j c )/<j c , A/x in accordance with 
Eq. @, we find 



4 > Albeit we review disordered state, we can consider that by parity of reasoning with the thermodynamical systems the 
microscopical susceptibility decreases from the maximum value x(°c) down to zero within the ordered region. Thus, the 
ergodicity region of ordered state is bound by the dotted line like the one shown in Fig. 3. 
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A = B(p)(a-a c ), 



B = 



1,1 c 1 f Xo 

1 + o CT XO - T — 

2 6 V<7o 



n 2 , Xo 

Xo og 



n -i 



(53) 



The coefficient B(p) diverges at the point corresponding to divergence of derivative d\/da at a = a c (see Fig. 2). At 
fixed value of thermal disorder it is necessary to make the expansion in — q n in equality (^0|). Then taking into 
account the dependence (kt9) in linear approximation, we obtain 



A = go c [(p/p c f 3 - 1 , 
where the critical value of memory parameter meets the ergodicity breaking point p c . 



(54) 



V. FLOW OF SAND GRAINS MOTION ON INCLINED SURFACE 



The viscous flow of sand grains on plain trajectory y — y{x) will be shown to represent the simplest example of 
self-organized criticality. The time dependencies of horizontal x — dx/dt and vertical y — dy / 'dt velocity components, 
and slope of sand surface y' — dy/dx fix the behavior of a system. In the autonomous mode they satisfy the equations 
of Debye relaxation 



dx 
dt 



x 

TV 



(55a) 



dy 
dt 



V 

jo) ' 



(55b) 



cV 
dt 



y'o-y' 

5 

ts 



(55c) 



where the t x , Ty°\ t$ are corresponding relaxation times. In Eqs. (55) is supposed that the stationary state meets 
the sand grains being at rest (x = y = 0), the slope y' ^ plays a role of control para met er. 

Since the motion along different directions is not independent we ought to add to Eq. ( |55a| ) the tran sverse component 
/ = y/j of force caused by motion along y axes (7 is the kinetic coefficient). As a result, Eq. (55a) takes the form 



T T X 



-x + a 1 y, 



(56) 



where a = "f/r x . Note, that in accordance with diffusion equation y — Dy" , where D is the diffusion coefficient, the 
stated force is proportional to curvature of sand surface: 



/ = {Dh)y". 



(57) 

In the stationary state, when x = 0, equation (|56| ) solution yields the dependence y = ax+const for the tangent 
to sand surface. On the other side, here the friction force ( |57| ) becomes proportional to longitudinal component: 
/ = t^x. 

Taking into consideration the relationship (|57]) and obvious equality dy 1 /dt = y' + y"x, the expression ( 55c ) assumes 
the form 



rsy' = [y' - y') - (ts/d) yx. 

Analogously for the vertical component we obtain 



y 1 • 

—yx, 



1 



1 

Jo) 



1 



/ (0) ' 

y T v 



(58) 



(59) 



Here the nonlinear addends of higher order are omitted and relaxation time t v is introduced that depends on stationary 
slope y' . 
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Equations (|56|). (|58|), and ( |59| ) coincide with the Lorentz system (4), if the horizontal x and vertical y velocity 
components respect to the order parameter r\ and conjugate field h, and the slope y' - to the control parameter S. 
Then the coupling constants of the system (4) take the form 

A v = aT 1 , A h = T y / Tx , A s = t s /D. (60) 

Taking into account the relationships A v oc vt x , Ah oc vt v , and As oc vts obtained from the comparison with the 
microscopic model ([[]), it is seen that the values of phenomenological parameters 7 oc D oc t x oc v^ 1 are fixed by 
microscopic one v. In the previous section the self-organized criticality has been shown to be realized, if v exceeds 
the value v c . It implies that avalanche formation on the sand surface takes place spontaneously, if shear viscosity 
being proportional to the relaxation time t X) is bounded from above by critical value, This condition is similar to 
transition criteria of viscous liquid into turbulent flow mode. 



VI. CONCLUSIONS 



As the above shows at fixed values of thermal and quenched disorder intensities a, p the behavior of self-organizing 
system is represented, on the one side, by order parameter field 77 and fluctuation amplitude tp of conjugate field, and 
on the other side - by the couple of Grassmannian conjugated fields ip, tp. Within the microscopic representation the 
Bozon and Fermion gases, interacting between themselves by means of potential v in Hamiltonian ([!]), meet them. 
Since with transition to self-consistent field scheme the couple of three-tail vertexes, every meeting the v, forms the 



four-tail one whose anharmonicity parame ter w is in bare superpotential fll8|), it is possible to suppose that the relation 
w = v 2 is fulfilled. Then the expression (|A.q) f or the critical value of thermal disorder S c , in which we should take 
into consideration the factor kj -1 (see after (|18|)), takes the form 

(T\ 2 (C 2 



2 t v 



Sc = ~ - r- , (6i) 



where we pass to dimension quantities, T is the noise intensity reduced to temperature for thermodynamic systems. 
The obtained S c (v) dependence has the same character as in the last formula (|^). Identifying they, w e fin d the 
expression for disagreement parameter e that provides the time irreversibility in energy balance equation ( |A.2| ) : 

tq\ 2 2 / a \ 2 t 2 _ 2tt?l 



r2 



U T h TS T 



ro = ^r- (62) 



Here the dimension units are used, h is the Planck constant. The coefficient c is defined by ratio of adiabaticity 
parameter (a/<^) 2 , expressed in terms of spatial scales, to the corresponding value t^ts/t^ in terms of characteristic 
times. It is possible to believe that c is the constant, and disagreement parameter ( J62| ) is determined by square of 
ratio of quantum fluctuation time tq to macroscopic time t v of the order parameter change. Obviously, the condition 
e <C 1 is satisfied always. 

The peculiarity of self-consistent supersymmetic scheme presented in sections [n]- |y] is that it allows us to go out 
of framework of the adiabatic approximation. On conditions of its applicability, v <^ v c , the critical value of thermal 
disorder is defined by equality (|6l|). With growth of interaction parameter v the complete suppression of disordered 
state takes place, when values v overcomes a critical value v c = 0.252. In this end, the effects of mutual influence of 
superfield components manifest themselves substantially, and the dependences type of represented in Figs. 1-6 would 
be used. 

We have shown for the simplest example of sand flow on the inclined surface (section [v]) that above scheme 
represents the avalanches formation in self-organized criticality phenomena Jl(|. If can suppose that the effective 
potential form as usually the complicated landscape in system's configuration space || similar to the spin glass [ p0[ . 
Therefore the theory of diffusion over nodes of hierarchical tree, formed by set of statistical ensembles of nonergodic 
system, would be used for the complete description of self-organized criticality [ 2l|,^2| . The above approach pretends 
to describe the conditions for single avalanche formation of lowest hierarchical level only. As is seen from Fig. 4 one 
can determine the two system types: at supercritical values of interaction parameter v > v c the ordering process 
realizes independently on external conditions; in opposite case v < v c the system passes into self-organization mode, 



Obviously, the relationship 7 oc D of kinetic and diffusion coefficients represents the Einstein expression. 
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if the thermal disorder intensity is above the critical value er c . Examination of the simplest sandpile model shows that 
on the macroscopic level the condition of spontaneous avalanches formation is similar to turbulence criteria. 

Phase diagram of a system defines the total landscape pattern and allows to predict the self-organized criticality 
character as a function of external conditions. So, Fig. 3 shows that system is nonergodic if the quenched disorder 
intensity exceeds the maximum value p m fixed by value of v. According to dependence represented in Fig. 5 at 
supercritical value v > v c the self-organized criticality process does not requires the quenched disorder at all. It is 
worthwhile to note that numerical experiments Ref. J^]- ||, where the avalanches intensities are stochastic values, 
meet the large values of quenched disorder. 



VII. APPENDIX 



To giv e the form of synergetic system (6) for the field equations (14) we multiply Eq. (14c) by ip to left-hand side, 
and Eq. (14d) by ip to right-hand side, and add the results. Then for quantities 



S = ipip, 



{Wip)ip - ipVip 



(A.l) 



we obtain the continuity equation S + V j = 0, that, obviously, expresses the entropy conservation law for conserved 
systems. Respectively, combinations (A.l) of Grassmannian fields determine the entropy S and its current j. It 



is characteristic that Eqs. (14c), (14d) for Grassmannian conjugated fields ip, ip differ only by sign in front of the 

time derivative, so that dependencies ip(t), ip(t) coincide at it inversion. Namely this circumstance p rovid es the 

entropy conservation condition albeit for each of the fields ip, ip this condition is not fulfilled: according to (14c) in the 

t 

homogeneous case the quantity ip(i) decreases exponentially with decrement < _1 J Vq'' ' {r){t'))dt' , and the conjugate field 

o 

ip(t) increases with the equal value of increment. For the entropy S = ipip the pointed out processes are compensated, 
and magnitude of S is conserved. As a result, the obtained continuity equation does not contain feedback with order 
parameter ry. 

To switch on this feedback we ought to take into account the disagreement of right-hand sides of Eqs. (14c), ( p_4d| ) 
that reflects the macroscopic time irreversibility. In this aim we introduce the coefficient 1 + e in the right-hand 
side of Eq. (14c) defined by the disagreement parameter e 1 (its value is determined in section VI). Then in the 
right-hand side of continuity equation the term — eVq" {rf)S appears. In addition we take into consideration that self- 
organization process realizes only at stationary current j which leads to thermostat entropy S increase with constant 
velocity — Vj = (r r) /rs)S° (at that the entropy of self-organizing system AS = S° — S decreases naturally). As a 
result, entropy balance equation reads 



5' 



1 /t s )S° - eVo"(v)S. 



(A.2) 



From here in the stationary regime S — 0, the equality S — (T n /eTs)S° /Vq' (r?) is obtained that is reduced to the 
form (0) for the bare potential 



Vo=r) 2 /2 + r) 4 /12. 



(A.3) 



Now let us examine equation (14b) for the fluctuation field ip(r, t). In contrast to conjugate field h the fluctuational 



one if is the nonhomogeneous even in the stationary state. We use the approximation V 2 p = (t;/a) 2 (p, where a is the 
scale of stationary fluctuation variati on, £ > a is the correlation length to be the scale of the coordinate r. Then in 
the stationary state tp — from Eq. ( |l4b| ) the relationship follows 



ip - 



v '"( v )s 



^/af-V "{ V y 



(A-4) 



showing that the values tp, S are connected at expanse of anharmonicity Vo (f)) ^ of bare potential only. In addition 
the stability condition <p > requires that the scale of nonhomogeneity does not exceed the value (Vq" (rj))^ 1 ^ 2 ^. At 
much rigorous requirement (et/£) 2 Vo"(?7) 1 the addend — Vo"(r]) in the denominator may be omitted. The pointed 
out inequality meets the adiabatic condition that, however, relates not the time but spatial scales. 

In terms of adiabatic approximation Eq. (14b) for amplitude (p(t) of most probable fluctuation takes the form 



ip = -(Z/a) 2 cp + V '"(r))S. 

Taking into consideration S(r]) dependence in the stationary state, ip = 0, we obtain 

ip = (a /2)Vb'" (r]) /Vq" (rj) , a = S°/S c , where the characteristic value of control parameter is introduced 



(A.5) 



14 



S c = (e/2)(Ts/T v m/aY 



(A.6) 



Using approximation ( A. 3 ) we see that found dependence tp{rf) is reduced to (Q). 

Lastly, let us consider equation (14a) for the order parameter field ?y(r,t). Inserting there the dependence p(r)), we 
arrive at the Ginzburg-Landau-Khalatnikov equation 



V ■ 



V 2 ?/ = -dV/d V , V(r}) = V (r]) - (cr/2)lnVb"(?7) 



(A.7) 



that differs from (N) by gradient term appearance. For dependence (A. 3) the synergetic potential assumes the form 



V{rj) = V{rj) + r/ 4 /12, V( V ) = (1/2) ( V 2 - a In (1 + V 2 )) 



(A. 



differing from (||) by the addend rf /\2. This distinction is caused by the circumstance that the potential V is defined at 
constant field h, whereas V - at fluctuation amplitude tp fixed. In other words, the first potential is the field h function, 
whereas the second one depends on fluctuation amplitude p. The quantities h, tp represent the couple of conjugated 
stationary state parameters (like the volume and pressure in thermodynamics) and the synergetic potentials V(h), 
V(tp) are coupled by Legendre transformation V — V ~ hip. The state equation governing the h(tp) dependence follows 
from the condition h = — dV /dp. But it is made simpler to introduce the field (j> — —Vq' (t)) + V 2 r)= — 8Vo{r)}/5r) 
that is reduced to the force / in Eq. (0) with accuracy to gradient addend. Then equation (14a) assumes the form 



fj = <f> + ip. 



Comparing Eqs. (A. 9) and (6a), we find the important relationship 

h = T) + (<p + p) = (r) - Vo'(ti) + V 2 ?y) + p. 



(A.9) 



(A.10) 



In the stationary state (?) = 0) the amplitude of the most probable fluctuation tp = — (f> coincides with generalized 
force with accuracy to sign, and field h — r\ is reduced to the order parameter. In the general case the disagreement 
(f> + p =/= results in variation of order parameter in the course of time and difference between the fields h, p is 
conditioned by nonlinear component of generalized force </>. 

Notice that obtained Eqs. (A.7), (A. 5), (|A.2| ) and the Lorentz equations (6) coincide in their mathematical structure 
only. So, the entropy balance equation ( |A.2|) contains the negative feedback similar to equation (|6c|), however it is 
expressed by addend — eVo (r])S that is proportional to th e en tropy, whereas the corresponding term —Agr/h does 
not contain entropy. The differen ce o f addends Vo'"(rj)S in flA.5| ) and A^iyS in (|6b|), me eting the positive feedback, is 
less essential (for bare potential (AJ3) they are equal at all). Finally the field equation (A.7) for the order parameter 
differs from the Landau-Khalatnikov equation (||) by account of spatial nonhomogeneity only. As was noted, the 
physical reason for above distinctions is that the Lorentz equations contains the field h conjugate to order parameter, 
whereas the amplitude of most probable fluctuation ip is in the initial field equations (14). Since the h and tp play 
a role of conjugated parameters of system's stationary state, then developed field formalism and Lorentz scheme are 
mutually supplementary approaches - the first one is used at fixed value of most probable fluctuation amplitude tp, 
and the second one - at fixed field h. Obviously, the second case is realized much naturally (to pass to this we should 
exploit the state equation ( A.10| )). 
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CAPTIONS 

to the paper "Field theory of avalanches formation" by A.I. Olemskoi and A.V. Khomenko 

Fig. 1. Dependence of microscopic memory parameter q on thermal disorder intensity a at anharmonicity parameter 
w = 0.01 (dashed line meets the unstable state, a c is the ordering point). 

Fig. 2. Dependence of microscopic \ an d macroscopic xo susceptibility on thermal disorder intensity a at unhar- 
monicity parameter w — 0.01 (dashed and thin solid lines meet the unstable state; thick line - stable state; a c is the 
point of ergodicity breaking, <r c is the ordering point). 

Fig. 3. Phase diagram: (a) at w = 0.01; (b) at w = 0.02 (0 is the ordered phase, D - the disordered phase; E - 
ergodic, N - nonergodic; the dashed line shows the precise boundary of ordered ergodic region). 

Fig. 4. Dependence of thermal disorder intensity a c corresponding to the ordering point on anharmonicity parameter 

w. 

Fig. 5. Dependence of maximal intensity p m of quenched disorder at the boundary of ergodic region on anhar- 
monicity parameter w. 

Fig. 6. Dependence of nonergodicity parameter A: (a) on thermal disorder intensity a (curves 1-6 meet the p = 
0.6; 0.8; 1.0; 1.4; 2.0; 4.0, respectively); (b) on quenched disorder intensity p (curves 1-5 meet a — 0.0; 0.1; 0.2; 0.3; 0.4, 
respectively). The anharmonicity parameter w = 0.01. 
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